Polaron formation as a genuine nonequilibrium phenomenon 
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Solitons and polarons at nonequilibrium steady states are investigated for the spinless Takayama 
Lin-Liu Maki model. Polarons are found to be possible only out of equilibrium. This polaron 
formation is a genuine nonequilibrium phenomenon, as there is a lower threshold current below 
which they cannot exist. It is considered to be an example of microscopic dissipative structure. 
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Nearly 30 years ago, the nonmetal- metal transition by 
doping halogens in polyacetylene was discovered by Shi- 
rakawa et.al. 0]. As a result of observing a mobile de- 
fect in polyacetylene, Su, Schrieffer and Heeger Q stud- 
ied solitons within a lattice model (SSH model). Then, 
Takayama, Lin-Liu, and Maki Q showed that its con- 
tinuum version, the so-called TLM model, allows one to 
study solitons analytically. Subsequently, polarons were 
numerically found in the SSH model and Brazovskii- 
KirovaH and Campbell-Bishop [|| have independently 
demonstrated the existence of polarons within the TLM 
model. Subsequently, its electronic and transport prop- 
erties have been extensively studied both experimentally 
and theoretically Q, and it is known that the charge 
carriers in conjugated polymers are localized excitations 
such as solitons and polarons. In spite of these devel- 
opments, there still remain several important issues to 
be settled. Recently, the critical dopant concentration 
for the nonmetal- metal transition was improved [H and 
the absence of the soliton contribution to the current 
was reported Q. An energetically most preferable state 
among those with solitons, polarons, bipolarons, and no 
localized excitations is not well understood E3 ■ The 
stability of polarons against various perturbations such 
as an external electric field [HI, EH, 0, EH an< ^ thermal 
noise 0, E?| is still intensively studied even now, as is 
the dependence of the polaron velocity on the applied 
field [II El Qj| and/or Coulomb interaction [2(j. 

As of yet, only the destructive role of electric field and 
temperature on the stability of polarons has been stud- 
ied. In this Rapid Communication, we report a construc- 
tive role of current within the TLM model. Namely, we 
show that, at nonequilibrium steady states (NESS), cur- 
rent induces polarons that are absent at equilibrium. It 
is a genuine nonequilibrium property, since there exists 
a lower threshold current below which the polaron can- 
not exist. Although only the TLM model is discussed 
here, our analysis covers a wider class of systems. This 
is because, in the mean-field approximation, the TLM 
model is equivalent to the XXZ model, the extended 
Hubbard model and the Gross-Neveu model in field the- 
ory [ULliil. 



The new polaron solutions are possible both for the 
spinful and spinless TLM models. Thus, to emphasize 
the very role of the current, we discuss the spinless TLM 
model which is known to admit no polarons at equilib- 
rium. The self-consistent conditions out of equilibrium 
are derived based on the scattering-theoretical charac- 
terization of NESS proposed by Ruelle[Il, [H, HI- 

The Hamiltonian H = Hs + V + Hb is composed of 
Hs for the finite TLM chain, Hb for the reservoirs, and 
V for their interaction, which are given by 
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where ^(x) = (d(x),e(x)) is the two-component spin- 
less fermionic field satisfying a boundary condition; 
d{0) = 0, e(£) = 0, A(x) is the lattice distortion, n(sc) 
is the momentum conjugate to A(x), a kv (y = L,R) 
are the annihilation operators for reservoir fermions with 
wave number fe, hw kv represents their energies measured 
from the zero-bias chemical potential at absolute zero 
temperature, o~ x and a y are the Pauli matrices, £ is the 
length of the system, v is the Fermi velocity, A is the 
dimensionless coupling constant, and coq is the phonon 
frequency. We assume that the coupling matrix elements 
v k as well as the density of states of the reservoirs are 
energy independent [27[; thus, the integral 



dk 



w — ^>kv — i0 



it / dk\v k \ 2 5(u-uj klJ ), (y = L,R) 



becomes a positive constant T. 

Next we describe the mean-field approximation. Since 
we are interested in NESS, the self-consistent condition 
is derived from the equation of motion for the lattice 
distortion, 
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Namely, the self-consistent equation is written 

A(x) +irhv\(¥(x,t)a x y(x,t)) A i F -0 , (2) 

where A(x) is the mean-field NESS average of A(x), and 
(''')oo F represents the mean-field NESS average. The 
mean-field NESS corresponds to the initial state where 
two reservoirs are in equilibrium with different chemical 
potentials, and it is characterized as a state satisfying 
Wick's theorem with respect to the incoming fields ctkv 
[y = L, R) of the mean-field Hamiltonian, and of having 
the two-point functions (28l. |29| : 

{ a \u a k'y)oo = f„(tiwkv)5(k - ft') , {v = L,R) 

where afc„ corresponds to the unperturbed field a^, 
f v {x) = l/(cxp{(x — [i u )/T} + 1) is the Fermi distribu- 
tion function with temperature T and, chemical potential 
Hl = —eV/2 and [Ir = eV/2 (the Boltzmann constant is 
set to be unity). 

At first, we briefly review the previous results on the 
uniformly dimerized case [2^], in which the average lat- 
tice distortion is constant: A(x) = Ao, the fermionic 
spectrum has a gap 2|Ao|, and Ao obeys the gap equa- 
tion 
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where lo c is the energy cut-off. This reduces to a well- 
known expression at equilibrium in the absence of a bias 
voltage. This equation is valid when the chain length £ 
is sufficiently long. The average lattice distortion A is 
found to be a multi-valued function of the bias voltage 
when T < T* ~ 0.5571 x T c . But, in terms of the current, 
which is given by 



Ao I < I el Khuj c 



de ^ e \ , A » [f R (e)-f L (e)] , (4) 



it is a single-valued function at every temperature. In 
the above, Go = e 2 vT / {nh(v 2 + T 2 )} is the conductance 
in the normal phase. Thus, the temperature and the 
current are chosen as control parameters. The phase di- 
agram on the J-T plane and the current dependence of 
the average lattice distortion are shown, respectively, in 
Fig. [T]and Fig. OH (left), for A" 1 = 2.4. In these figures, 
the average lattice distortion, the temperature, and the 
current are scaled, respectively, by the zero-bias lattice 
distortion A c = ftcj c /coshA _1 at T = 0, the zero-bias 
critical temperature T c = 2fi,w c exp(7— \^ 1 )/tt (7: Eulcr 
constant), and the critical current J c = GqV c at T = 0, 
where V c = 2hio c exp(— A _1 )/e is the critical bias voltage 
at T = 0. The multi- valued property of the average lat- 
tice distortion with respect to the voltage results in NDC 
for T <T*. We note that NDC was reported in the field- 
driven SSH model [1^, as well as in the stochastically- 
driven XXZ and extended Hubbard models [U, [13]. How- 
ever, we also note that the NDC reported in these mod- 
els is different from the one discussed here. Since, as 



detailed above, the current is a multi-valued function of 
the chemical potential difference. This point merits fur- 
ther investigation. 

Next we investigate the solitons and polarons. Observ- 
ing that the only difference between the NESS and equi- 
librium cases is that the Fermi distribution is replaced 
with the averaged distribution {/i(e) + fn(e)}/2, the 
self-consistent Eq. ^) is expected to have similar solu- 
tions to the case at equilibrium. It is easy to verify that 
Eq.([2]) admits a soliton solution [33| similar to that of the 
equilibrium cased, 

A(x) = Ao tanh k s (x — a), k s = Aq/(Hv), 

where the amplitude Ao is the solution of the gap equa- 
tion ([3]) for the uniformly dimerized phase, and a = 0(£) 
represents the center of the soliton. At the same time, a 
midgap state appears with energy hiu = in the fermionic 
spectrum. Note that, even when solitons exist, the cur- 
rent is still given by (£f]). Then, following Brazovskii- 
KirovaQ and Campbell-Bishop Q, we look for a static 
polaron solution of the following form 



A(x) = Ao — hvnQ(t + — t-) 



U 



tanh kq (x cl i xq ) , tanh2ftoXo 
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where Ao and xo are parameters that are determined self- 
consistently, and a is the position of the polaron center on 
the order of I. As in the equilibrium case, the correspond- 
ing fermionic spe ctrum consist s of continuum states with 
energy hui = zty 7 (hvk) 2 + A§ (|fc| < uj c /v), and midgap 
states with energies hut = ±^/Aq — (Hvkq) 2 = zLHuib- 
Even though the coupling between the midgap states and 
the reservoirs is exponentially small for long chain length 
£, it still controls the occupation of the midgap states at 
NESS[34[. Therefore, one should carefully take a long 
chain limit, resulting in a self-consistent equation l[2"jl 
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where (3 = 1/T, Is is a contribution from the continuum 
states , and Is is a contribution from the midgap states 
with energy \Huj\ < |A | (see [Hj]). Comparing term by 
term, the gap equation ^ is obtained, and the equation 
for energies ±Hlub of the midgap states 



lob dio 



sinh H/3lu 



A \/h \Juj 2 — A 2 /h 2 (uj 2 — lo 2 b ) coshhflu + cosh 



sinh h[3uj b 



2vk cosh H/3u>b + cosh ^ 



(5) 



3 




FIG. 1: Phase diagram on the J-T plane. 
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FIG. 2: The left figure shows current- voltage characteristics 
at T = 0.05T C . The right figure shows the current dependence 
of |Ao| at T = 0.05T C . In these figures, the solid line is stable, 
and the dashed line is stable only at a constant current. In 
the right figure, only the bold solid line admits polarons. 



Note that the current is still given by (H} for the po- 
laron solutions. Eqs.® and §5§ have a nontrivial solu- 
tion only when the current (equivalcntly, the bias volt- 
age) lies between the lower and upper threshold values 
Ji(T) < J < J 2 (T) (Vi(r) < V < V 2 (T)), and the tem- 
perature is lower than T* , under which the system shows 
NDC. As seen in the left figure of Fig. [3l the polaron 
width 2xq and amplitude A = 2(Hvko) 2 / (\Aq\ + Hujb) 
are decreasing functions of the current. When the cur- 
rent (equivalcntly, the bias voltage) approaches the lower 
threshold Ji(T) (Vi(T)), the polaron width diverges and 
the polaron amplitude approaches the soliton amplitude 
2|Ao|. This indicates that the polaron splits into a 
soliton-antisoliton pair. On the other hand, when the 
current (the bias voltage) approaches the upper thrcsh- 
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FIG. 3: The left figure shows the current dependence of the 
amplitude A (the solid lines) and the soliton size 2xo (the 
dashed lines) at T = 0.05 x T c (the thin lines) and T = 0.2 x T c 
(the bold lines). 2xo is scaled by L = hv/ A c . The right figure 
shows a typical lattice profile at T = 0.05 x T c . From top to 
bottom, J= 10~ 3 J C , 10~ 5 J C , 10" 10 J C and 10" 15 J C 
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FIG. 4: Current dependence of the positive-bound-state en- 
ergy hu B at T = 0.05 x T c (the solid line) and T = 0.2 x T c 
(the dashed line). 



old J 2 (T) (V 2 {T)) : both the width and amplitude of the 
polaron vanish, and the polaron solution reduces to the 
uniform solution. Typical profiles of the polaron solution 
are shown in the right figure of Fig. [3] 

As mentioned above, |Ao| is a multi- valued function of 
the bias voltage and, for a given voltage, several uniform 
phases are possible. Although this suggests the possibil- 
ity that collective local excitations can separate uniform 
domains with different values of |Ao|, there exist only 
those interpolating uniform phases with the same |Ao|, 
such as the solitons and polarons just discussed. This 
is because charge conservation implies that the current 
J remains constant over the chain, and Ao is a single- 
valued function of J. Also, it is interesting to note that 
the existence of the polaron solution is related to the lin- 
ear stability studied previously p9j. Indeed, the polaron 
solution exists when the uniform phase with Ao is stable 
both at constant current and constant bias voltage (the 
solid curves in the right figure of Fig [2J , but it does not 
exist if the uniform phase is unstable at constant voltage 
(the dashed curve in the right figure of Fig [2]). Because 
of this property, there is one-to-one correspondence be- 
tween the current and bias voltage intervals where the 
polaron solution is possible, Ji(T) < J < J 2 {T) and 
Vi(T) < V < V 2 (T), respectively. This aspect and the 
non-existence of the polaron solution for T > T* deserve 
further investigation. Note that the states on the thin 
solid curve do not admit polaron solutions. 

The possibility of the polaron solution at NESS can 
be qualitatively understood as follows. Recall that the 
polaron at equilibrium is possible only in the spinful 
case. With the corresponding fermionic state, the lower 
midgap state is occupied by two fermions with oppo- 
site spins, and the upper midgap state is occupied by 
an unpaired fermion. In the half-filled spinlcss case at 
equilibrium, such an asymmetric occupation is not pos- 
sible. This seems to suggest the necessity of the particle- 
hole symmetry breaking for the polaron formation. This 
seems to suggest that it is necessary for the particle- 
hole symmetry to break for polaron formation. In con- 
trast, at NESS, the particle-hole symmetry is broken 
by the bias voltage even for the half-filled spinlcss case. 
This is because the fermionic occupation is controlled by 
(/i( e ) + /fl(e))/2, which is not symmetric under the ex- 
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change of particles and holes. 

It is interesting to note that, at low temperatures, 
the width J%(T) — Ji{T) of the current interval that 
admits polaron solutions increases with an increase in 
temperature, while the width V2 (T) — V\ (T) of the volt- 
age interval decreases with temperature. These be- 
haviors of the current and voltage are consistent, be- 
cause the phases admitting polaron solutions tend to 
become insulating phases as T — > 0, which implies 
lim T _ +0 (J 2 (T) - J 1 (T))/(F 2 (T) - Vi(T)) = 0; thus, the 
decrease of V 2 {T) — V\(T) with an increase of T does 
not contradict the increase of J2(T) — Ji{T). Because 
of the discontinuity at T = of the R.H.S. of Eq. (©, 
which behaves like the Fermi distribution function, ab- 
solute zero temperature is a singular point. Indeed, at 
T = 0, Eq. and Eq. (J5J) admit a polaron solution with 
hjj B = (tt 2 /16 + 1)~ 1/2 A C only when 
F = {(7r 2 /16 + l)cosh 2 A- 1 }- 1 / 2 exp(A- 1 )^ c and J = 0. 

The existence of solitons and polarons has been verified 
by spectroscopic experiments, where the energies of the 
associated midgap states arc observed 0, HjJ |3g|. Fig. d] 
shows the current dependence of the energy Hub for the 
midgap state at T = 0.05 x T c , 0.2 x T c (< T*). As shown 
in the figure, Hub is a monotonically increasing function 
of the current, and it approaches for J — > Ji(T); and 
Ao for J — > J2 (T) ; this reflects the change of the polaron 
profile. Polarons in a spinful system possess this same 
feature, since the corresponding self-consistent equation 
is obtained simply by replacing A in ^ with 2A. Namely, 



the energies ±Hujb of the midgap states associated with 
NESS polarons change from to ±|A | as the current in- 
creases, while those with equilibrium polarons in a spinful 
system are fixed at ±1tlub = ±|Ao|/-\/2- Such a current- 
induced shift of energy spectra might be observed by 
spectroscopic experiments. 

In summary, we have studied solitons and polarons in 
the open spinless TLM model, and in particular, we have 
shown that polarons are possible only out of equilibrium. 
The polaron formation is a genuine noncquilibrium phe- 
nomenon, as there exists a lower critical current J\{T) 
(equivalently, a lower critical bias voltage Vi(T)), below 
which polarons are not possible. This observation sug- 
gests that the new polaron is an example of microscopic 
dissipative structure. Also, we have shown that the 
critical temperatures for polaron formation and the ap- 
pearance of the negative differential conductivity (NDC) 
agree, although polarons are not allowed at the current 
found in the NDC regime. The energies of the midgap 
states associated with polarons are shown to crucially de- 
pend on the current, which might be observed by spec- 
troscopic experiments. 
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